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ABSTRACT 


By adapting the streamline model developed by Caudle and LeBlanc 
for flow through porous media, a streamline model for in-situ combustion 
was developed. 

Initial tests were made to compare the streamline model for in-situ 
combustion with actual field trials conducted by oil companies. The results 
were extremely close. Asa result, in-situ combustion experiments may 
be examined and optimized to obtain the best possible results prior to im- 
plementation in the field. 

Conducting numerous tests by changing assumed values indicated 
that areas in the process must be re-examined in order to determine the 
economic feasibility of using in-situ combustion to recover the residual 
oil remaining in the nation's reservoirs, kerogen from oil shale, and the 


hydrocarbons from tar sands. 


lv 





TABLE OF CONTENTS 


PAGE 

N OWLE DCEME NTS аа ааа... n iii 
КОК АО ЛТ ООО een een iv 
ME OF CONTE NO ......................................... V 
DE TABLES. A l.cs vi 
ШИШ FIGURES AA ss ео, ера ная vii 
CHAPTER 

] ШЧТЕО ШЕ ОПМ......................................„... 1 

II ТАЕ ЛИЕ ШЕ МОРЕ... ааа... 3 

III ТІМ-БПИФФКТВПТЕӘТІОУМ МОПЕТ,......................... 9 

IV EXPERIMENTAL INVESTIGATION ....................... 26 

V CONCLUSIONS AND RECOMMENDATIONS ............... 54 
APPENDICES 

A FEON DO cr r nenene aana 97 

B ООМРОТЕНК РЕВОСКАМ.................................. 62 

C КО КШ ОСТАТОВЕ........................................ 68 
RA o nad 74 
VITA 





LIST OF TABLES 


TABLE PAGE 

I Recovery Data for Basic Streamline Model of 

In-Situ Combustion in an Unbounded System ............... 30 
H Recovery Data for Streamline Model Utilizing 

Moss Relationship Between Advancement of 

Combustion Front and Injected Air FluX................... o 
III Recovery Data for Basic Streamline Model of 

In-Situ Combustion in a Bounded System .................. 46 
IV Characteristics of Sun Oil Company Field Test............ 48 
V Characteristics of Magnolia Petroleum Company 

PIC. ШМ .-..........4......›.. 50 


vi 





LIST OF FIGURES 


FIGURE PAGE 
1 Unbound Five Spot Pattern Showing Generated 
streamlines from Wells Having Equal Rates............... 7 
2 Typical Streamlines ina Two Well System................ 8 
3 ро o A сое 8 
4 Gate's View of the In-Situ Combustion Process............ 10 
5 Allred's View of the In-Situ Combustion Process .......... 11 
6 In-Situ Combustion as Used in the Modified 
ре ООЛО ......................................... 12 
1 Effect of Injected Air Flux on the Combustion 
ОИТ асе ERST иеа ааа 14 
8 Process Zones in a Моп-Ііпеаг буѕіегт.................... 21 
9 Process Zones in a 1.1пеаг бузїет......................... 21 
10 Streamlines in the Basic Reservoir Used for In- 
ЕЕ О аео ны о a re Sv 28 
11 Location of Oil Zone in Basic Test at Breakthrough ....... 31 
12 Combustion Front Location Using Martin and Moss 
Relationship Between Injected Air Flux and Com- 
ри ЕП ОЛІ Каб(е оі Айуапсв............................. 92 
13 Effect of Variation in Hydrocarbon Consumption........... 35 
14 Effect of Variation of Zone One Mobility................... 37 
15 Effect on Oil Production by Variation in Injected 
М о ee O 39 
16 Effect on Time of Breakthrough by Variation in 
Da О UN a ee 40 
17 Combustion Front Location at Breakthrough with 
iamiam nsiDprected Air BlIUxC er e eee mme 42 


vii 





LIST OF FIGURES (Continued) 


FIGURE PAGE 

18 Typical Pattern Showing Imaging Technique, ............... 43 
19 Location of Combustion Front in a Bounded 

System o o ае e 222272222 45 
20 Streamline Model Burned Zone versus Sun Oil 

AA ЕЕ... ао 49 
21 Streamline Model Burned Zone versus Magnolia 

Petroleum Company Field TeSl............. ne 91 


vili 





CHAPTER I 
INTRODUCTION 


In 1972 the United States became aware of an energy crisis. The 
crisis took the form of shortages of fuel for production of electricity and 
heating fuels for homes and businesses. The causes are many and varied, 
but may be traced to one reason — economies. Since the start of oil pro- 
duction, the cost per barrel has been kept so low by artificial barriers that 
the return of the dollar invested in oil or gas exploration and production is 
below the return of investments in other areas. Thus the energy crisis is 
not the lack of fuel per se, but the lack of economic fuels. 

All persons intimately related with the oil industry know that existing 
means of production may optimistically recover fifty percent of the oil in 
the nation's reservoirs. Thus one could conservatively estimate that over 
90 billion barrels of oil still remain in established reservoirs, since the 
total production of the United States from 1920 to 1969 was just over the 90 
billion barrel mark (1). In addition, potentially 1,781 billion barrels of 
kerogen from oil shale and 399 billion barrels of oil from tar sands exist in 
the United States (2). Based on an estimated rate of consumption of 46 mil- 
lion barrels per day by 1980, the above resources could supply the United 
States for over 135 years. 

The fallacy with the above conclusion is that the price of a barrel 
of oil would have to rise drastically to make production of the residual oil, 
oil shale, and tar sands profitable if present production concepts are main- 


tained. For example, present methods of oil production from oil shale is 





by mining the oil shale and retorting the shale above ground. Not only is 
the process environmentally distasteful, but the resulting cost per barrel 
of oil is between four and five dollars (2). 

During the past ten to fifteen years, one process has received exten- 
Sive investigation —— in-situ combustion, often called fire flooding or under- 
ground retorting. Through in-situ combustion, residual oil in existing oil 
fields could be recovered. Synthetic oil from oil shale and tar sands could 
be recovered without strip mining, and possibly at lower costs than the 
present. 

To the author's knowledge no attempt has been made to simulate in- 
situ combustion by use of the Caudle-LeBlanc Streamline Model (3). Ex- 
pressed simply, the Caudle-LeBlanc model is a mathematical representation 
of the movement of one typical particle of oil. The following study describes 
the basis of the model and the modification of the model to include the specific 
characteristics of in-situ combustion, with the concluding sections comparing 
the model results with the results of actual field tests conducted by various 


organizations. 





CHAPTER II 
THE STREAMLINE MODEL 


The theory and development of the streamline model has been re- 
viewed by many different authors (3, 4). For this reason only a brief review 
of the model will be shown. Following this brief introduction the model 


modifications will be shown. 


The Basic Streamline Model 
The basic equations were developed from Darcy's Law and included 


the following assumptions: 


negligible gravitational effects inside the reservoir 
incompressible homogeneous fluids 


horizontal system, infinite in extent 


BD фо N m 


homogeneous, isotropic media 


The most general form of Darcy's Law is as follows: 
q = > me —_—_——— (2. 1) 


]t is also known that the fluid flux is equal to the flow rate divided by the 


area through which it flows, thus Equation 2.1 may be written as 


E E o (2. 2) 





4 
In order to determine the streamline path it is necessary to know its velocity 
inthe x, y, and z directions. Since it is known that the velocity is the fluid 
flux divided by the porosity of the material through which the fluid is flowing, 


Equation 2.2 may be written as: 


Velocity = v = =- . k P (2.3) 
ф Фи de 
The second important equation is the continuity equation. It is a partial 
differential equation derived from the material balance around a point in 


three-dimensional space. This equation is written as follows: 


a (pux)  a(puy)  3(pu) ao 
OX 78 ч 07. E^ dt. BIS 


When the volumetric flux components of Darcy's Law, Equation 2.2 in the 
x,y,z directions, are substituted in Equation 2.4 the following equation 


results: 


2 pk, 9Ф з [оку 99 4 [0k, 99 : 0 a5 
— |————— г ---- -------- + ——— = — 
OX | u 0X oy | и ду д> \ и 02 ot (2.5) 


Applying the assumptions made earlier, Equation 2.5 becomes: 





2 2 

др др 
=e = 0 (2.6) 
OX ду? 


which has become known as LaPlace's Equation. By transforming to cylin- 
drical coordinates and allowing radial flow around a point, the equation may 


be solved for the flow potential in a horizontal, isotropic, homogeneous 
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porous medium. Restoring the resulting equation to rectangular coordinates 


and allowing for a multiwell system, Equation 2.6 becomes: 


n 
p(x,y) = Py - a ia n| (x-x;)? * G-y9? | (2.7) 


where Pms the mean reservoir pressure. 
By differentiating Equation 2.7 with respect to x and substituting the 
resulting derivative into Equation 2.3 the velocity at any point in a rectan- 


gular system becomes 





- q (x-xi) 
-Xj 
"o Ca Жр E 
TON i21(x-x)^ (у-у) 
In a similar manner, the velocity inthe y direction becomes 
n 
1 а (у-у;) 
vy(xy) = - ) (2. 9) 


In order to follow the particle along a particular streamline, the dis- 
tance of travel must be known. Given a finite amount of time, the distance 


a particle travels at a constant velocity may readily be calculated as 


ds = vd or AS = VAt (2. 10) 


Assuming a particle was initially located at point (xj, yj), the location of 


the same particle after time At could be calculated as 





(2011) 
V NNI vy ^t 


By starting on the wellbore of the injection well and systematically 
reiterating Equations 2. 11, a trace of any particle's movement may be made. 
Figure 1 shows the paths of particles, from a single five-spot pattern with a 
central injection well. The time to reach a specific point may be determined 
by the summation of the At's required to cause sufficient movement of the 
particle to reach that specific point. 

In a given two-dimensional fluid system, an infinite number of par- 
ticle paths (streamlines) could be used to describe the flow from a source to 
a sink. However, to simplify computation, a representative number of 
streamlines may be chosen to represent all fluids emanating from the source. 
А few of the streamlines from a source are shown on Figure 2 (5). By 
dividing the distance between two adjacent streamlines into two equal parts, 
the streamline may be said to represent the flow of all particles within the 


dividing lines. A streamtube has thus been defined and is shown on Figure 3. 











generated streamlines 


ее” 
^ 
e — 


FIGURE 2 Typical Streamlines in a Two-Well System 


flowtube 





generated streamline 


FIGURE 3 Typical Flow Tube 





CHAPTER III 
IN-SITU COMBUSTION MODEL 


Prior to developing the modifications necessary to convert the stream- 
line model to a streamline model for in-situ combustion, a discussion of the 
process of in-situ combustion is warranted. 

In-situ combustion is not a new process. Recovery by a combustion 
or heat wave process was patented in 1923. F. A. Howard was granted a 
patent on a process in which air and a combustible gas were pumped into an 
injection well and ignited (6). The method involves ignition of the formation 
in an injection well, followed by propagation of a combustion front through 
the reservoir. Combustion is maintained by the injection of an oxygen- 
containing gas, such as air, to react with reservoir hydrocarbons. Asthe 
flame front, which is really glowing embers of the reservoir material, pro- 
gresses through the reservoir, oil and formation water are vaporized, driven 
forward in the gaseous phase and recondensed in the cooler part of the for- 
mation. In turnthe condensed fluids push oil into the producing wells (7). 

Since 1923 field experiments have been conducted by the Russians in 
1935 (8) and others in Oklahoma (9) and Kansas (7), with the most recent 
test occurring in the Green River Oil Shale deposits by the U.S. Bureau of 
Mines (11, 12). The objective of the early tests and those occurring in the 
later years was to recover the residual oil left after primary drive and 
secondary methods had been utilized. In the 1950's in-situ combustion ex- 
periments were conducted on oil shales with the intent of retorting the 


shale underground as an alternative to mining the shale and retorting above 
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13 
ground. Other tests were conducted in highly viscous oils with the intention 
of reducing the viscosity of the oil in order to make the oil more mobile and 
thus producible in a smaller period of time (10). 

The complex nature of the in-situ combustion process, with it's 
many inter-related displacement mechanisms, makes exact interpretation 
of field data impossible. This is understandable considering that even two- 
phase flow of fluids through a porous homogeneous medium, which has been 
studied for many years, cannot be adequately described. Asa result, as- 
sumptions must be made based on experience and the results of many trials. 
Gates and Ramey (10) conclude that, as a result of their interpretation of 
field tests in California, Figure 4 represents the in-situ combustion process. 
Allred's (13) beliefs are somewhat different and describe the in-situ com- 
bustion process as shown on Figure 5. With the exception of the gases 
shown in Figure 5, both representations are the same. This author has 
simplified the two figures to enable utilization in the streamline model simu- 
lation. The in-situ combustion process as utilized in the simulation is 
shown on Figure 6. 

Along with the composition of the various zones, the velocity of the 
fronts is important. Since the burn front is the driving mechanism in the 
process, its velocity is the basis for determining the velocities of the other 
zones. Two experiments indicate that the injected air or gas must reach a 
certain rate before movement of the front will begin. Figure 7 represents 
the results of the aforementioned experiments (9, 14). Further experiments 
may show that the rate of advancement of the combustion front depends on 
the oil content of the reservoir, reservoir porosity or some other variables. 
What is important at this time is that the relation between air flux across 


the burn and the rate of movement of the burn does exist. 
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Other such relations and assumptions had to be established in the 
modification of the streamline model, but to ensure their complete compre- 
hension they will be discussed as they affect the streamline model modifi- 
cation. 

In developing the original streamline model we assumedthat the fluid 
was incompressible. The effect of thisassumption is extremely close to reality. 
Merrick (18) compared the results of flow of real gases, ideal gases, and 
incompressible fluids through a porous media and determined that the dif- 
ferences in results were extremely small. This assumption will also be 
made for the in- situ combustion model of flow in a reservoir, but to utilize 
the relationship of the advancement of the combustion front and the injected 
air flux, the injected air flux must be in surface volumes. To properly 
change the injection rate to surface volumes we must look at the gas laws. 


For an ideal gas 


pV - nRT (3. 1) 


where V isthe volume of a given number of gas moles n at absolute pres- 
sure p andtemperature T.  R isthe universal gas constant. For a non- 


ideal gas the relation is 


pV = ЕТ (3.2) 


where z isthe compressibility factor ofthe gas. The compressibility is 
a function of the pressure inthe reservoir pq. Thus, Equation 2.3 de- 


scribing an ideal gas becomes 


V = = 2m (353) 
$ Pm 
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where Zn is the compressibility factor of the fluid injected into the reser- 
voir, and p,, is the average reservoir pressure in atmospheres (standard 
pressure is one atmosphere) and the surface and reservoir temperature are 
equal. 

Equations 2.7, 2.8, and 2.9 must also be modified since they 
were developed for noncompressible flux. With an ideal fluid, the volume 
at standard conditions equals the volume at reservoir conditions. Since we 


are using a compressible fluid, the quantity at surface conditions becomes 


Ise zm Ч тез (3.4) 


Since it is also known that the mobility of a fluid is equal to its permeability 
divided by its viscosity, Equations 2.7, 2.8, and 2.9 become 


n 


Z 2 2 
р = Pm” EE Lasci Фл (ве + (у-у;) | (3:5) 
а (х-хі) 
Zm 1 
ee) SSC =e om ee. (3.6) 
“Sx 2т\фрү, à і (хех)? + (yy) 
n 
zn (у-у) 
= yg Lisco, 7 в . 
"Sy 2rhópq, D n (х-х;)< + (y-y;) (3.7) 


From the above equations the pressure and velocity of the standard fluid may 
be determined. The standard fluid in this thesis will be the air injected into 
the reservoir. It is called the standard fluid to keep the process as general 


as possible. In addition, the standard fluid has particular significance in the 
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conductivity ratio which will be explained in the coming pages. 

To determine the amount of production, the velocity of the burn is 
necessary Since, as will be shown later, the movement of all the zones 
containing oil are dependent on the movement of the burn. Equations 3.6 
and 3.7 will generate the velocity of the standard fluid, which in this case 
will be the air injected into the reservoir. Knowing the velocity of the air, 
and by the relationship of the air velocity and burn front velocity as expres- 
sed in Figure 7, the velocity of the burn may be computed. From Figure 7 


the velocity of the burn is 


Eb - mus. + С (3.8) 


W = SLv. +C (3.9) 
y y 

where SL is the slope of the line and C is the line constant. Note that 
at approximately thirty feet per hour, the relation between the burn and the 
standard fluid changes. At this point the values of SL and C will be 
changed. In the streamline computer program these values are used as input 
data to provide as much flexibility as possible. As previously mentioned, 
insufficient testing has been completed to prepare exact figures for the burn 
front movement with regard to the rate of injection of air, but the fact that 
the relation does exist must be considered. 

The initial streamline model developed in Chapter II was for a single 
fluid. For a multifluid condition, such as a waterflood or in-situ combustion 
where different fluids exist, the concept of the conductivity ratio must be 
applied. Caudle (15) presents a detailed explanation of the concept and thus 


it will not be discussed here. The conductivity ratio is determined as fol- 
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lows: 


р; -Р 


eL (3. 10) 
L (Pg-1 Z Pg?) 
Na y [E 
ІІ Àk 


k=1 


where y is the conductivity ratio, p, is the pressure at the source or in- 


jection well, p, is the pressure at the sink or production well, A, is the 


p 
mobility of the standard fluid, and à is the mobility of the fluid in the k 
zone which has a pressure p, at its interface with the zone preceding the 
k zone. L isthe number of different zones. 


Substituting Equation 3.3 in Equations 3.8 and 3.9 and multiplying by 


the conductivity ratio 


Pm 
ЖШ. № ae Y) C (3. 11) 
X X Zm 
Үр, = SL (vs i 2 + С (3. 12) 
y У Zm 


As a result the velocities (x and y directions) may be computed for 
the burn front. The velocities of the other three fronts will be similar, but 
will differ by some factor due to the accumulation of oil in front of the burn. 
By referring to Figure 8 it may be seen that the velocity across point B 
will yield a quantity Qp of a known amount if the cross-sectional area of 
point B is known. By the same rationale a fluid moving across point 2 
will yield a quantity of fluid Qo of known amount if the area and velocity 


of point 2 are known. Additionally, if a constant rate of injection is estab- 
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lished, the quantity of fluid passing point B will, in some way, be related 


to the quantity of fluid passing point 2. Mathematically, 


Әр = Vb А, (3. 13) 
Qo a Ao (3. 14) 
Qo = Е20, (3415) 


Thus, by substituting Equation 3. 18 in Equation 3. 15 


Considering the constant rate of injection, the following conclusions may be 


made: 


Q. = vs, Ay (3.17) 


©. Vso Ao (3. 18) 


Equating Equations 3.17 and 3. 18 and solving for Ao | 


А 


"Sp b 


(3. 19) 


А 
2 
ус 


2 


Substituting Equation 3.19 in Equation 3. 16, 


Vs, Ab 
Vo = ы = Роу, Ар (3. 20) 
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Solving for the velocity at point 2, 








= Е P 
Va 5 Еҙур - (3.21) 
Sp 
By similar reasoning point 3 on Figure 8 would have a velocity of 
Vs 
3 
Vg = Fav E (3.22) 


As shown on Figure 6, the fourth zone is composed entirely of gas, and by 
a piston-like displacement the zone is pushed by the air injected into the 
reservoir. Thus the velocity of zone 4 is that of the standard fluid at zone 
4 as computed by Equations 3.6 and 3.7. 

Considering all the x and y components of the velocities of zones 


2,3 the movement of these two points may be computed as follows: 


в V 
270,75 
X 
bx 
Еур У 
2 by Soy 
Vo = E (3.24) 
y 
"by 
Еу; V 
2S D.'s 
уз = E 5 ox (3.25) 


V 
"bx 
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FIGURE 8 Process Zones in a Non-Linear System 
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FIGURE 9 Process Zones in a Linear System 
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F 
у “Зу 
уд = ee (3.26) 


To compute the values of Fo and Ға , a movement of the burn must 
be considered. Referring to Figure 9, a movement of the burn from B to 
B' will cause the displacement of all oil existing in the area between B and 


B' that is not consumed as fuel in the combustion process. Thus 


OC 
Oil moved = Vp (^r - 6) (1.0 - ы". һы 


(1.0 - Sg.) 
(322) 
uf d Vp (Gp - $4) (1.0 = fp) (1.0 - Sg, - Sy.) m 
en (1.0 - Sg.) 
Since 
v; = A (3. 29) 
Ai Sort 
d Vg (dp - 6) (1.0 - ,) (1.0 - Sg; - S) ЕД 
er eee ; 
і dt A óc. (1-5) 
Since 
dV, /A; 
‘py. = 3Yp/ ^ (3.31) 


1 dt 
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Equation 3.30 may be written as 


(65 -06.)(1.0- £ )(1.0- S5, - S,,.) 
у. = ЕСІ ee (3.32) 
бег (1.0 - Sg; | 


where Фе is the effective porosity through which the fluid must move and 


is computed as 


бен = 9, (1.0-S,-S,) (3.33) 


As was discussed earlier, the velocity of the burn is proportional to the 
velocity of the air flux across the burn front. Thus if we consider point 2 


as the i-th point, 


V Ус 
pe 2 (3.34) 
us "s, 
Ys, Vb 
“bo = = (3 39) 
Sh 
and Equation 3.32 becomes 
(05 7 6, (1.0 - ) (1- 8g, - Sy) so Vb 
== - (3.36) 





V9 = A < ЕН 
фен (1 - Sgo) Y 


Comparing Equation 3.36 with 3.21, we findthat F, is 
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(фр - 6,) (1.0- ao 
Е, = ы лаш” (3.37) 


9, 1 - Sg.) 


By mutual displacement, the movement of point 2 to 2' has moved 
point 3 to 3' by an amount equal to F. Thus the movement of point 3 to 
Br is 

V3! "b 


Уз! = Fo E" — (3.38) 
sb 


The movement of point 3’ to 3'' caused by the movement of the burn and 


not the movement of point 2 to 2' is given by Equation 3.32 as 


(or - 5) (1.0 - $) (1.0 - S - буу Y Vs, Y 


53! 3! S3! 





uer = 
V 
A T Sp 


(95029) 


Combining Equation 3.38 and 3.39 to obtain the total movement of point 3 





WS", 
(ép - 65) (1.0 $) (1.0 - 8g, - Syy,) V, Vp 
dert) Sp 
(3.40) 


Comparing Equation 3.40 to Equation 3. 22, Е. is found to be 
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- 1203 : 2 
(6p - $9 (01.0 - 5) (1-8, - S) 
Ез = > _ _ _ t F 
Фен (1-56) 


(3.41) 


Substituting the equations developed in this chapter into the streamline 
model developed by Caudle and LeBlanc (3) will result in an in-situ combus- 


tion streamline model. 





CHAPTER IV 
EXPERIMENTAL INVESTIGATION 


Sensitivity to Assumptions 

As was discussed earlier, little is known about the in-situ combustion 
process. Present technology prevents us from obtaining cores from the 
reservoir as the process moves through it. Asa result, numerous assump- 
tions had to be made. The first assumption is the relationship between the 
advancement of the combustion front and the injected air flux. The effect on 
production and time of breakthrough by changes in this relationship must be 
determined. The second assumption is the amount of reservoir hydrocarbon 
consumed in the propagation of the combustion front. As may be seen in 
Equations 3.36 and 4.40, the velocity of zones 2 and 3 are directly re- 
lated to the amount of reservoir hydrocarbon consumed. The third variable 
is the mobility of the various zones within the process. 

In order to determine the sensitivity of the streamline model to vari- 
ations in these assumptions, a basic test was conducted with the following 


data. 


1. The rate of advancement of the combustion front in relation to 
air flux is in accordance with the curve in Figure 7 established by Martin 
(14) because it shows a realistic slowing of the combustion front as the air 
flux declines. 


2. The saturation of gas, oil, and water in each zone is 
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Zone 1: So 1700 
Zone 2: Sg EUN 
So = ЗА 
E = qu 
Zone 3: So = 0.50 
So е 0.30 
E 20220 
Zone 4: Sg = 1,00 


3. The rate of injection was established at 600 scf/hr. 

4. The mobility of each zone was assumed to be 25, 11, 12, 25, 25 
darcy/centipoise for Zones 1, 2, 3, 4, and5 respectively. The wisdom of 
using these mobility ratios in light of the saturation and make-up of each 
zone may be dubious; however, tests were conducted to determine the effect 
of changes in zone one mobility only. The others were not examined due to 
reasons to be explained later. 

9. The amount of oil consumed in the propagation of the combustion 
front was assumed to be ten percent of the total hydrocarbons in the reser- 
voir. 

6. The average amount of oil in the reservoir was assumed to be 
27 gallons of oil per ton of reservoir material. 

7. Pattern utilized was an unbounded inverted five spot, 100 feet 


square. 


The above data were incorporated into the computer program shown 


in Appendix B. The plot of the streamlines is shown in Figure 10, and the 
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production data are shown in Table I with the location of the oil bank at 
time of breakthrough in Figure 11. In subsequent tests the location of the 
combustion front will be compared with the basic model plot in order to 
show the amount of the reservoir swept by the burn. 

The relationship between the advancement of the combustion front 
and the air flux is fundamental to successful simulation of in-situ combus- 
tion. The two curves showing this relationship in Figure 7 have been pub- 
lished as a result of experimental studies (9, 14). Asa result, these curves 
must be assumed to include such effects as oxidation potential of the oil, 
heat conduction and heat convection. The basic program was computed 
using the relation developed by Martin (14). To determine the significance 
of the difference in the Martin and Moss (9) curves, the second test was to 
use all the data from the basic test with the Moss curve. Figure 12 com- 
pares the results of the tests, where the basic data using the Martin curve 
is a solid line and the second test using the Moss curve is a broken line. 
From the start of ignition to 129 hours, the fronts are propagated radially. 
The positions of the fronts shown in the upper right quadrant is the position 
of the combustion fronts at breakthrough. The number of hours shown in 
that quadrant is the time of breakthrough of the indicated front. Table I 
shows the production data for the second test. The difference in the break- 
though time and the total production may be explained by review of Figure 7. 
At high air flux rates, the front of the basic test moves at a faster rate. 
Production data confirms this in that the initial oil production is 97 percent 
greater in the basic test, indicating a larger area swept in the same amount 
of time. As the air flux rate drops due to the expanding surface area of the 
burn front, the Martin curve declines more rapidly than the Moss curve. 


After 629 hours the combustion front of the basic test has reached the pro- 





TABLE I 


RECOVERY PERFORMANCE FOR BASIC TEST 


Gas Produced Oil Produced Water Produced 
(ft?) (bbls) (bbls) 
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Basic Test (Martin Relationship) 
Second Test (Moss Relationship) 
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FIGURE 12 Combustion Front Location Using Martin 


and Moss Relationship Between Injected 
Air Flux and Rate of Frontal Advance (14, 9) 





TABLE II 


RECOVERY DATA FOR SECOND TEST 


Gas Produced Oil Produced Water Produced 
(ft?) (bbls) (bbls) 
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duction well and the streamlines midway between two production wells have 
stopped. Ihe reason for these streamlines stopping may be found in Equa- 
tion 3.10 which has been expanded below to include all zones considered in 


this thesis. 











Pi 7 Bp 
y = 
ND Pp Ho O PADS 
e + — + + -————- 
SI Ag ^2 dg 4 А5 


(4,1) 


As the streamlines between two production wells move out from the injection 
well, the pressure at those points decreases which in turn decreases the 
conductivity ratio. Since the conductivity ratio is directly proportional to 
the velocity of the burn, as shown in Equation 3. 11, a decrease in the con- 
ductivity ratio will cause the combustion front to slow and eventually stop. 
The same would be true for the streamlines reaching the production wells, 
except that, as each zone is produced, its term in the denominator of Equa- 
tion 4.1 is dropped thus causing an increase in the conductivity ratio. The 
slowing and stopping of the streamlines between the production wells causes 
the oil front ahead of the burn to stop before reaching the production wells. 
The next area of investigation was to determine the effect of varia- 
tions in the amount of hydrocarbon consumed in the propagation of the com- 
bustion front. Equations 3.38 and 3.42 show the relationship between the 
velocity of the different fronts to the amount of reservoir hydrocarbon con- 
sumed. As the amount of reservoir hydrocarbons consumed is increased, 


the velocity of each front will be decreased. Figure 13 shows that as the 
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amount of hydrocarbon consumed is increased from 10 percent to 60 percent 
the amount of production and the time required for breakthrough increase. 
When the consumption rate is set at 70 percent of the reservoir hydrocarbon, 
the total production and time of breakthrough decrease. 

The results indicate that, as the velocity of the combustion front de- 
creases, the amount of production also decreases. Confirmation of this 
effect should be obtained when the injection rates are varied. However, 
prior to studying the effect of variation in injection rates, the last assump- 
tion should be discussed. 

The third assumption made in the utilization of the streamline model 
is the mobility of zone 1. As the combustion front moves across the reser- 
voir, the water and oil have been driven forward or been used as fuel for 
the combustion process. In shale the clay would shrink much as a brick does 
during firing. Logically, the permeability of the zone behind the combustion 
front should increase and thereby enable the injected gas to pass through the 


zone more easily. The mobility of this zone is 


ми (4. 2) 


where k isthe permeability of the material through which a gas or liquid of 
a specific viscosity (u) is passing. With an increase in permeability, the 
mobility should also increase, assuming the viscosity of the gas or liquid 
does not change. The mobility of the first zone has the greatest effect on the 
velocity of the combustion front since it is the first zone which contains the 
standard fluid. Equation 3.10 shows the rationale for the above statement. 
Figure 14 shows the results for four trials using a mobility of the first zone 


of 25, 50, 75, 100 and 200. As may be seen, the mobility of zone one has 
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little effect on the time to breakthrough, but does affect the amount of oil 
produced. For the streamline approaching the production well fastest, the 
conductivity ratio changes as the most advanced zone is produced. Equation 
4.1 was expanded to include all zones considered in this paper. As may be 
seen, if the leading zone is produced, the conductivity ratio expression de- 
creases by one term in the denominator since it no longer exists in the 
reservoir for that particular streamline. Since the mobility of zone one is 
the same as the standard, the conductivity ratio approaches unity as each 
zone is produced, regardless of the value used as the mobility of zone one. 
Thus the mobility of zone 1 has little effect on the streamlines that travel 
the shortest distance to the production wells. The effect on the production 
may also be explained by Equation 4.1 in the same way as the stopping of 
the burn front was determined. However, the larger the mobility of the first 
zone, the larger the value of the denominator and the smaller the conductivity 
ratio. Thus, the greater the value of the mobility of zone 1 the shorter the 
distance from the injection well to the final position of the burn front. 

The mobility of the other zones may also change during the in-situ 
combustion process, but their effect should not be so significant as a change 
in zone 1 mobility since they are filled with a greater percentage of liquid 


and the temperature decreases rapidly between combustion front and zone 2. 


Injection Rate Variation 

Additional tests were conducted to determine the effect of increased 
injection rates on the total production and the time of breakthrough of the 
combustion front. For these tests, the relationship between injected air 
flux and the advancement of the combustion front was as developed by Martin. 


The results of these tests are shown on Figures 15 and 16. These results 
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FIGURE 15 Effect on Oil Production by Variation in 
Injected Air Flux 
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FIGURE 16 Effect on Time of Breakthrough by Variation 
of Injected Air Flux 
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were rather interesting in that at higher injection rates the total production 
was lower. The explanation of this may be found in Figure 7. At high in- 
jected air flux rates the combustion front moves rapidly between the source 
and sink. The movement is much faster than the case of lower injection 
rates and results in an early breakthrough. As a result, the streamlines 
off the direct line between the injection and production wells are not able to 
move very far from the injection well. Thus, a smaller amount of oil is 
produced. The opposite action occurs at lower injection rates. The stream- 
lines on the direct line between the injection well and production well move 
rapidly, but do not breakthrough at the production well before the other 
streamlines have moved a considerable distance from the injection well. 

The graphic evidence of this action is shown in Figure 17 where the broken 
line represents the location of the combustion front for an injected air flux 
of 4000 scf/hr, andthe solid line represents the combustion front for an 
injected air flux of 600 scf/hr at breakthrough. 

The results of the tests on the effect of increased injection rates in- 
dicate that economic considerations must be made to determine the maximum 


injection rate versus the time and amount of oil recovered. 


Bounded System 
Heretofore, the reservoir has been considered infinite in length and 
width. Now the reservoir will be examined with a boundary existing just 
outside the production wells. The pattern boundary is shown in Figure 18. 
In previous tests we noted the decrease in oil production due to the 
fingering effect and rapid expansion of the combustion front, thereby forcing 
oil out of the pattern area thus making it non-recoverable. As is shown on 


Figure 11, a majority of the oil within the pattern has been pushed away 
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--- Injection Rate of 4,000 ft? /hr 


— Injection Rate of 600 {t>/hr 





FIGURE 17 Combustion Front Location at Breakthrough 
with Variation in Injected Air Flux 
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from the injection well and in fact away from the production wells. In order 
to close the system, the imaging technique will be utilized. The imaging 
technique is a well-known method for confining streamlines inside a de- 
sired boundary. The technique is described in detail by Collins (13). In 
cases of constant reservoir thickness only the flow rate and location of each 
well is necessary. 

The pattern used, including the image wells is shown on Figure 18. 
The location of the burn front at breakthrough is shown on Figure 19. Com- 
paring the bounded system which has the same properties, injection and 
production rates as the basic model shows that the start of production and 
breakthrough occur fastest in the bounded system, and that the production 
rate of the bounded system lags that of the unbounded system. This could 
be expected since the imaging technique prevents flow across the established 
boundaries which are shown by the broken line in Figure 18. As a result, 
the oil laden zones are directed towards the production wells; however, in 
the unbounded system, the zones ahead of the burn take a more circuitous 
route as shown in Figure 1 and sweep part of the reservoir outside the area 
of the test. 

Reviewing the production data in Table III for the bounded system 
shows the start of production and breakthrough of the combustion front occur 
earlier than in the basic unbounded system. The bounded reservoir pro- 
duces about one-third as much oil as the unbounded. The importance of the 
existence of natural boundaries is therefore apparent. The existence of 
faults close to the pattern could cause similar reductions in production as 
the artificial boundary utilized in the above test. Fortunately, Lin (16) has 
developed a technique for construction of boundaries for the streamline 


model that may be used to simulate the shape of natural barriers found in 








TABLE III 


RECOVERY PERFORMANCE FOR BOUNDED SYSTEM 


Gas Produced Oil Produced Water Produced 
(ft) (bbls) (bbls) 
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reservoirs. By including boundaries, the area swept and total production 
may be determined for various locations of the injection and production 


wells. 


Comparison with Field Tests 

To establish the ability of the streamline model to determine the 
location of the combustion front, data from two field trials were used as 
input data for the model. The first test utilized information provided by 
the Sun Oil Company (17). The test was conducted in the May-Libby reser- 
voir which has the characteristics as shown in Table IV . The pattern used 
was an inverted five spot with dimensions as shown on Figure 22. Ofthe 
information provided on Table IV, only the porosity, injection rate, and well 
pattern and spacing were utilized. The rationale for using only this infor- 
mation will be explained at the end of this section. The saturation and mo- 
bility of each zone was assumed to be the same as that of the basic test in 
the first section of this chapter. In addition the amount of reservoir hydro- 
carbon existing and consumed in the propagation of the combustion front was 
the same as in the basic test. The burned area is shown as the shaded 
area in Figure 22. As may be seen, the area burned in the field test and 
that indicated by the streamline model are extremely close. 

The second test was a comparison of the streamline model and a 
field test conducted by the Magnolia Petroleum Company (9). Table V 
shows the characteristics of the reservoir and Figure 23 shows the dimen- 
sions of the inverted five spot pattern. As inthe previous test, the satur- 
ation and mobility of each zone were as used in the basic test as well as 
the amount of reservoir hydrocarbon existing and consumed. The burned 


zones are quite similar except in the region between the injection well and 





TABLE IV 


CHARACTERISTICS OF SUN OIL 
COMPANY TEST 


* Well Pattern: Inverted Five Spot 
* Well Spacing: As Shown on Figure 10 
* Average Porosity: 31. 200 

Depth of Pay: 3, 400 ft 

Thickness of Pay: 4.4 ft 


Average Reservoir 1,069 md 
Permeability: 


Average Amount of 740 stb/acre-ft 
Oil in Place: 


* Average Injection 8,000 cu. ft/hr 
Rate: 





ж 
Used as input for streamline model 
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— — — Sun Oil Company test in May-Libby reservoir 


“cer. Streamline model 


FIGURE 20 Streamline Model Burned Zone versus 
Sun Oil Company Field Test (9) 





TABLE V 


CHARACTERISTICS OF MAGNOLIA PETROLEUM 
COMPANY FIELD TEST 








* Well Pattern: Inverted Five Spot 


* Well Spacing: As Shown on Figure 23 







* Average Porosity: 27.2% 









Depth of Pay: DOS 





Thickness of Pay: 17 ft 


Average Reservoir 7,680 md 
Permeability: 


Average Amount of 1, 350 bbls/acre-ft 
Oil in Place: 


* Average Injection 3,200 scf/hr 
Rate: 


* 
Used as input for streamline model 
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== =» Magnolia Petroleum Co. Test 


Tee Sireamline Model 


FIGURE 21 Streamline Model Burned Zone versus 
Magnolia Petroleum Company Field 
Test (17) 
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production wells land 4. The authors of the article published on the burn 
indicated that the reservoir properties dictated the non-radial formation of 
the burned zone. The streamline model of in-situ combustion is based on 
a uniform permeability; therefore, if the permeability of the areas in 
question are lower than the rest of the reservoir, little of the injected air 
would go in the direction of wells number land 4. In addition, if the reser- 
voir was tighter between the injection well and production wells 1 and 4 
the production of these two wells would be lower than that of the other wells. 
The authors confirmed this fact. 

A third test was attempted using the information published by the 
U.S. Bureau of Mines with regard to an in-situ combustion experiment in 
oil shale (11, 12). The result was not satisfactory because of the nonuniform- 
ity of the permeability of the oil shale due to extensive fracturing around the 
injection and production wells. The test pattern consisted of an inverted 
five spot with the first set of four production wells approximately twenty- 
three feet away from the injection well. A second set of production wells 
approximately fifty feet away from the injection well were located on the 
Same axis as the first set. With the extensive fracturing consisting of 
electrolinking, explosive fracturing and by hydraulic fracturing, areas of 
highly permeable oil shale existed and enabled excellent communications 
between the injection and production wells. Since the permeability in the 
rest of the reservoir was extremely low, all injected air would flow in the 
direction of the least resistance, i.e., toward the production wells. Making 
the assumption that the oil shale was a homogeneous, isotropic deposit, as 
was an assumption for development of the streamline model, was therefore 
erroneous. Asa result, the test by the U.S. Bureau of Mines could not be 


Simulated by the streamline model. 
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The articles relating the two field trials did not contain all the in- 
formation necessary to reproduce the field tests. As was stated above, the 
porosity, well spacing, and injection rate were the only data taken directly 
from the field tests. This was done for numerous reasons. First, the 
streamline model assumed that the reservoir was only one foot thick. Thus, 
the combustion front could be assumed to be a vertical front. Although this 
assumption is valid for a thin reservoir section, Moss has found that the 
assumption cannot be made for a thick reservoir (9). Secondly, the perme- 
ability is tied directly to the mobility of each zone. Since this area is one 
which is unknown at the present time, the permeability of the reservoir is 
ofnouse. Lastly, the difference in production, as demonstrated by the 
comparison of the Moss and Martin relationships between rate of advance 
of the combustion front and the injected air flux, would prevent meaningful 
comparison since it is not known if the aforementioned relationship is the 


same for all reservoir materials. 


Discussion 

The comparison of the streamline model to field trials demonstrated 
the model, as it exists, can approximate the location of the combustion 
front. Some ofthe variables that affect the combustion front location have 
been exam ined and the effect determined. The effect of other variables that 
warrant additional study are listed in the next chapter. 

As mentioned earlier, no attempt was made to determine the loca- 
tion of the oil bank or to duplicate production data due to the complexity of 


variables necessary to determine the velocities of zones 2 and 3. 





CHAPTER V 
CONCLUSIONS AND RECOMMENDATIONS 


Conclusions 
Ás a result of this study, the following conclusions may be 


drawn: 


1. The streamline model can be used to study the effects of vari- 
ables on the process of in-situ combustion. 

2. The streamline model can be used to determine the location of 
the burnfront with an accuracy dependent upon the accuracy of the input 


data and the homogenity of the reservoir. 


Recommendations 


1. The relationship between the rate of advance of the combustion 
front and the injected air flux should be determined for the specific reser- 
voir materials in order to increase the accuracy of the streamline model 
for in-situ combustion in the specific reservoir. 

2. A study of the problem of determining the saturation of various 
zones is necessary to obtain more realistic values of the zone composition 
and its mobility. 

3. Larger pattern areas should be utilized to determine the best 
economic relationship between injection rate and pattern area. 

4. An expanded study to include variation in reservoir thickness 


would be of great value. 
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DO 


9. The effect of dip should be considered as it affects the flow of 
oil during the combustion process. 

6. The results of this study were based on a reservoir with uniform 
permeability. With the extensive use of fracturing, a study comparing the 


effects of permeability changes is warranted. 
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APPENDIX A 


FLOW CHART 


READ: Number of Streamlines; Well Radius; 
Reservoir Porosity, Pressure, Density; 
Oil Content, Thickness, Boundary; 
Well Locations and Rates; Mobility, 
Saturation of each Zone; Flux versus 
Advancement of Burn Data. 


Print Data 


Convert English Units 
to Metric Units 


Compute Pressure at 
Injection and Production 
Wells 
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Plot Location of Wells 













Locate Starting Point 
of Each Streamline 


Iterate For Each Injection Well 
Iterate For Each Streamline 


Locate Starting Point 
of Each Zone 


Iterate For Each Injection Well 


Iterate For Each Streamline 






DO 


DO 


99 


Combustion Front at Yes 


Boundary of Reservoir 


Compute Pressure and 
Velocity of Standard Fluid 


Determine Zone at 
Production Well 





Й 
Й 
fi 
| 


Compute Conductivity Ratio 


Compute Air Flux at 


Combustion Front 





Air Flux Below Minimum | Yes 


Compute Velocity of 
Combustion Front 





= 


1 


bo 


* 


Compute Velocity of 
Zones 2,3, and 4 


Yes 


Combustion Front at 


Production Well 


Compute Increment of Movement 





Compute New Location 
of Each Zone 


Determine Zone at 
Production Well 


Redefine Production Ratios 


Compute Incremental Gas, 
Oil, and Water Production 
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Plot Location of 
Combustion Front 


Store Zone Locations 





Iterate For Each Injection Well 


Injected Air Produced 





Compute Time of Movement | 


Test Each Streamline 


Compute and Print Recovery 


Performance 





END 
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3 
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APPENDIX B 


COMPUTER PROGRAM 


PROGRAM COMR(TNPUTsOUTPUT) 

THIS PROGRAM PLOTS THE LOCATION OF THE BURN FRONT 

DIMENSION @ (7) Хы (7) ¢¥w7) МОЕ, (41 } 9 Т50М (7941) 6х (7 +41) беа 7 , 
со (7+5) ) «Си (741) 9С (7041) СҮ (7041) *XEB(UT4 51) » YEB(TS^4]) v 
TBl (CT $51) sXF3(Q4 1) » YFB(A41) V ATPROUT $4 1) 9 ХЕ2 (7:41) з ҮЕг (7, 
41) 5 ХЕе (41) э\Р2 (4 } ) ЭСАЗРЕО (7,41) 9ХС3 (7,41) 9УЕЗ (7,41) 9 
ХЕЗ (41) 9\ЕЗ (4] ) эхЕ4 (7.4] ) 4, YES(7$ 413 e XF 64 (41) «YF 4 (41) 

PEAD 3309 NST RAPINIS NDATIRI FMB PURI ZM PH 

READ 3909 PVOLMX9Hs XMAX 9 YMAX 9 IPY 

READ ^30, SLOPEOs*BONE*FLUXM 

HEAD ^40, SLOPET*ETWOsF LUXMN 

READ 580, HHO*FUELs*GALPF TN 

READ 520, RATIOZ+RATIOS9RATIOG/WRÁTIZ9WRATIZIWRATIG+GRATIZIGRATIA 
GRATI4 

ЕО. ООО 

PRINT 500» NST,NP.SNISNUATSRISPCD4,£M,PM 

PRINT 495, РУО МХ,» Н, ХМАХ» УМАХ, Іру 

PRINT 450» SLOPEOsBONE 

PRINT 530» SLOPLT»BTWO 

PRINT 420» RHO+FUEL+GALPTN 

PRINT 470» RATT[02«$RATIO3*RATIOA 

READ 540, SMB,TMB,FORMB,F IVMB 

PRINT 3509 FMB+SMB+TMB)FORMB+F ]VMB 

1 TOTzNP«NI 

READ 3455 (QUI) sX& (Iyo Yu (ID) , Iz], IT OT) 

READ 5450, ӘМІМ,ОМІМКР,АКЕА»5рТ5Т 

PRINT 570 

PRINT 380» (0 (1) Хи (Т) зүм (1) 612], 1ТОТ) 

PRINT 460» GMIN*GMINPs:^REA$SDTST 

SVALPV=AREAZH“ POR/SQMIN® 3600, 

SVALUEZIPV*SVALPV/NOAT 

CORRECTION FACTORS FROM ENGLISH TO METRIC UNITS 

GASPRÜ(JsK)z0.0 

WwATIPRD(J4,K)z0.0 

MDSTx=MpPSTAY30,48 

MOSTY-MCSTY*30,^G 

XMAXzXMAX* 30.48 

YMAXzYMAX* 04.48 

SBIST 07 1730,48 

Н=Н 30.48 

RI=SOTST/50.0 

NTOT=NI+NP 

PQ-55.720,48 

OMINZUMIN®7,88055 

OMINP=OQMINP#7.88055 

NIT ITSA\P 4] 

CTR=1.0 

PI-0,0 

PP=0,0 

NOIA i=j ITOT 

CORRECTION FACTORS FROM ENGLIS TO METRIC UNITS 

Хм (1) Ху (1) 30,469 

үм (1) =үр (12430,66 

0(1)-0(1)%/.88055 
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30 
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PISPI=(0(1)*ZM)/(1205664*FMB+H9pmM) “ALOG((XAW(NI])+RI=Y9(]1)) 42. 
) (УМ (м1) Ум (1) ) 92) 
РРЕРР- (0(1)*2М)/(]}2.5664%7МВ%НеРрМ) $ АЕОС ( (Хм (} ) « АТ-Хи (ТГ) ) “92+ 
1 (УИ (1) -Уи (1) ) %+2) 
CONTINUE 
CALL HGNPLT 
CALL PLT (10.015.593) 
CALL PLT (*3.0,=3.043) 
CALL PLT (7303.012) 
CALL PLT (3,01 3.012) 
CALL PLT (3.0»-3.0%2) 
CALUT PET (3:03: 012) 
SCL=.0010024788 
DO 25 Iz]1*1T70T 
CALL SYMBOL (Xw(I)*SCLsYW(I)*SCUL,*2191*0«2*71) 
CONTINUE 
DO 40 J=NIleNTOT 
HSL (CJ) Z&ST*Q(CJ) /uMIN«04. 45 
PRINT 370 
IBENSL (J) 
00 30 к-1.1В 
X(JsK)zRJ*COS(1*6,283c9K/NSL(UJ) ) «Aw (J) 
Y(JsK) zRIJ"SIN(1]*6,2832*K/NSL (CJ) b. eYw (CJ) 
XE2( JK) =X ТК) 
NECK USK yay (ein) 
ХЕЗ (К) =Х (К) 
УЕЗ (К) =\ (УК) 
XE& (CJ«K) XUJSK) 
УЕ 4 (К) =У (К) 
Ow (J K)=0eV 
COlJ4K)20.U 
С#(,),К)=0,0 
ТВТ (К) =0°0 
TSU =9. 0 
CONTINUE 
CONTINUE 
CONTINUE 
DO 290 J=zNI1»NTOT 
IRENSL (Y) 
00 299 Kz1]*1B 
IF (GW(JsK),GT.0,.0) GO TO 260 
ХЕВЕХ ЛК) 
YFE=SY (J.k) 
XFe-XE2(JsSK) 
YF2=YE2(J9N) 
ХЕЗ=ХЕЗ( уе К) 
YF 3zYE3(J*K) 
XAF4zXC4 (J34K) 
YF65zYE4( Ja) 
СЕРТКЕ К) СЕ, УК) + 5СьЗ) 
TSUMISTSUMIJAK) 
IF (^pnSiXFB) GT. XMAX) 270960 
CONTINUE 
IF (ABS(YFB) eGTeYMAX) 270970 
CONTINUE 
VELXP 70,0 
VELYR=0.0 
VELX2=0.0 
VELY2=0,0 
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VELX370,) 
VEUY 3-00 
VELX4=0,0 
VELY4=0, 
PR=0,U 
P2=0,0 
P3-0,0 
Р4=0.0 
PKSTzZM/(14., 5664 9FMB&8H*pM) 
VKSTzZM/(6«282329h«PORU PM) 
ПО 80 |=]з[Т01 
FACéG=(xp4-XAW(L)) 229 (YFG4-YW(L))+e2 
PL=P4=PKSTQ(L)*ALOGIFACI) 
VELY4AsVELY^«((VKSTeQUL) /CTR)* (YF 47 YW(L) ) /FAC4) /POR 
VELXaz*VELX4* CUCUKST*Q(L)/CTR) * (XF47XW (L2) ) /FACA4) SPOR 
VELS4-SGRTUVELXa * s 2e VEL Y 992) 
FAC3s(xF3eXU(L) ) *22« (YF 37Yw (L) ) «a2 
PZPI-PKSTV(L)*ALOGIFECI) 
VELX3-VrzLXJ« CCVKST*"Q(I(L) /CTR) * (XF 37XW (CL) ) /FAC3) 
VEEN G-vVonro* ( (vino) /CTR)@ (YE 3"YW(L))/FAC3) 
VELS2=SGRTIVELX3#%24VELYSt #2) 
FAC22(xrF2-Xw(L)) v*«2«(YF2-7 Yu (UL) ) e«2 
Pe-zP2-PKST*OUL) *XALOG (FACZ) 
VELX2zVELXC* (CCuyKST*Q(U)/CTRY)S (XF2*XW (L) ) /F ACO) 
vELYZ*VELYG*CUVRST9"Q(I(L) /CTR)* (YF27* Yw (LD) ) /F AC2) 
VELS?-SCRTÜVELXx2**32«yELYz Y «2) 
FAC8zUXFB9XWw(L))*e2e(Yr dm Yw(L) sec 
EPHePRePKST*OUD)vAEOOG (FACB) 
VELXPzVELXP« CCVKST9Q(L) /CTR) * (XFB7XW (L2) /F ACB) 
VELYP-VELYU* CCVKST*G(L) /CTRO) * (YrFB" Yw (CL) ) /F ACH) 
BU CONTINUE 
VELSRzSGRT (VELXxB* ?2«VyELYBS* 92) 
IF (SQRT (CABS (CXF 4) -PQ) **2* (&RSUYFA) -PQ) 2*2) BLT. RI) GO TO 90 
СТУ= (РР-РГ} И (Емб® ( ((28-РГ)ИЕМВ) + ( (Р2-РВ) ИЗМЬ) + ( (РЗ+Р2) ИТМВ) + { (Р4- 
] P3)/FORMB) + ((PP>P4)/FIVMB))) 
G0 TO 120 
90 CONTINUE 
IF (SQRT (CARS (XF3) P0) **2e (ABS (YF3) РО) * 92) JLT.RI) GO TO 100 
CTY=(PP=p1)/(FrB* (((p5>PI)/FMB) + ((P2-PB)/SYB) > ( (РЗ-Р2) ТУВ) < { (рр= 
1 P3)/FORMB))) 
Со то 120 
100 CONTINUE 
IF (SoRI((LABS(XF2) -PO) *$2« (ARSIYF2) -PG) 992) JLT,RI) СО ТО 110 
CTYz(PP-pI)/(FMB*(((pB8sPI)/FMB)«((P2-PB)/SMB)*«((P392P2)/TM*B) )) 
Go To 120 
110 CONTINUE 
DEMZFMBRB9(((PP-pB) /FMB)e((PB-PI)/SMB)) 
CTY=(PPp=-PpI)/DEM 
120 CONTINUE 
VELSA-VELSB9*PM/ZM*CTY 
IF (VELSA,LTeFLUXM) 1405130 
130 CONTINUE 
VELBZSLOPEO* VEL SA*BONE 
60 ТО 150 
140 CONTINUE 
IF (SLOPET*VELSA,LE.ABS(BTWwC)) GO TO 270 
.  VELBZSLOPET* VELSA*BTWO 
150 CONTINUE 
AMTOILzGALPTN*RHO*,.00006€684 
POR=ZPUR=AMIOIL 


0 






160 


180 


190 


200 


210 


бо 


FACTe-AMTOIJUL* (1, 0-FUEL)* VOLO2/(V0OLO2*V0OLO3) )/ (POR? vnLO2) 
FACT3-AMTOIL* (1,0 -FUEL) * (VOLO3/ (VULO2*VOLO3))/ (POR9 VOLO23) *FACT2 
РОВ-,25 
МЕГЕДСЕУЕГО/УЕТ 5В 
SVELX4=VELAG 
SVELYS=VELYS 
VELT4=VELS4 
IF (SORT (CABS (XF 3) -PQ) #24 (ARS (YF 3) -PQ)%#2).,.LT.RI) GO TO 160 
SVELX3ZVEL^3*VELFAC9FACT3/POR 
SVELY3-zVELYS3*VELF ACSF ACT 3/POR 
VELT3-SGRTISVELX3* 92e 5 VELY3«992) 
CONT PRUE 
IF (SaoRI (CCcABS(OXF2) PG) 92 (ARS IYF2) -PG) 8*2) ,LT.RI) GO TO 179 
SVELX2zVELA2*VELFAC"FACTZ/POR 
SVELY2zVELY2* VELF ACSFACTZ/POR 
VELTz-SonnRT(ISVELA29*2« SVELY2 3 82) 
CONTINUE 
SVELXB=VELARSVELFACIPOR 
SVELYPzVELYH* VELF AC/pOR 
VELTR-*SGRTISVELAB?* *2«SVELYB* 92) 
VELT!'zvEL T^ 
IF (УЕ Тао оТ, ЕСТМ) VELTU=VELTA 
IF (VELT3. GT, VELTM) VELTM=VELT2 
IF (VELT2eOTeVELIM) vELTr"zvELT2 
IF (VELTBeOTeVELIM) VELTV=VELTR 
IF (PB+6T.0,0) GO TO 180 
DO )80 Lz1sHP 
ВАО =59РТ ( (АВЕ ХИ (1) ) +6024 (УЕВ- уу (0) 0) 82) 
IF (FADeGTeRI) GO TO 180 
TBT CJ* K) ZTSUMI 
Qw(J4 K)zGMIN/NSL CJ) 
Go TO 2/0 
CONTINUE 
DT-RI/VELTM 
IF (TSUMI*DT. GT. SVALUE) 190,209 
CONTINUE 
DT-SVALUE- TSUMI 
CONTINUE 
XFB8zxFBRepT*SVELX8 
YFBzYFRepDT*SVEL YG 
XF2=xF2+*pT*SVELX2 
YF2=YF2+pT*SVELY2 
XF3zXF3«0T* SVELX3 
YF3zYF34pT"* SVELYÀ 
XF4zxF4«pT* SVELX4 
YFAzYF4«pT* SVELYA 
CALL PLT (XFB*SCLesYFBPSCL,2) 
YSUM]=TSUMÍ+DT 
5DTST=70+7*30+48 | 
IF (SORT (CABS (XFS) 2PQ) *92«* (ABSUYF A^) "PGQ)* 2G) JLT.RID) 60 ТО 210 
RATIC=0,0 
WRATJO=9.0 
GRATIO=0,0 
GO TO 2^0 
CONTINUE 
WRATIO=WRATI4 
RATIO-RATIO^ 
GRATIO=GRATI4 
IF (SQRFC(CARSUXF3) -PQ) * *2«* (ARS UYF3) -PQ) * *2) ,LT.RI) GO TO 220 
GO TO 2*0 





230 


240 


250 
260 


210 


280 
290 


300 


310 


320 
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CONTINUE 
KATIC=RATIVZ 
WwRATIOzwRATI3 
GRAT IOzGRAT I3 
IF (5СНІ ((АВ5(ХҒ2) “Р0)9%%2% (Ар5 (ҮҒ2) “РС)992),1 Т.НІ) GO TO 230 
G0 TO 240 
CONTINUE 
RATIO-RATIO2 
"RAT TIOzwRATI?2 
GRAT{O=GRATI2 
ADDPRO=CT* (GMIN* CZM/PM)/(NSLUJ) «VELSB*POR))*(,00000181* (VELSRB*pOR& 
CTYS 1PM/7M))+,00U399) 
CO(l J+K)=CO(J+K) + ADOPROY“RATIO 
GASPRD(Y,K) =GASPRD (Ugh) ¢ADDPROSGPATIO 
WATPRD( YG eKI EWATPRD (09K) tADDPROAHRATIO 
IF (ABSUXFB) , CT. XMAX) 2705259 
CONT INUE 
IF (ABSUYFB) ,GTe YMAX) 27053260 
CONTINUE 
rU TA NAE) 60 TO 70 
CALL SYvpgOL (XFB*SCL,YF B8" SCLS,07,11»04. »-1) 
CONTINUE 
AlJıK)=XFB 
У (У;к)ЕУЕВ 
XE2( YK) =XF 2 
YE2( UK) =Yre2 
ХЕЗ ЈК) = ХЕЗ 
ҮЕЗ( ЈК) = ҮРЗ 
XE4 CJsK)zXF 4 
YE (Js K) zYF 6 
TSUM(J, K) Z T SUM] 
CY(J, K) zCTY 
CONTINUE 
CONTINUE 
СОР= 0,0 
CGP=0+0 
CWTRP=0,0 
CTR=0+0 
CwWP=0 ® 0 
DO 3C0 JzNI]SNTOT 
IB=NSL (J) 
00 300 К-1%18 
CWP=CWPeQw (J9K) = (SVALUESTBT (Sok) ) %20000062898 
COP=COP+C0 (J9K)*.20000062898 
CGP=CGP+GASPRD (J9K) * „009935314 
CWTRP=CWTRP+WATPKDI(J» A 8 
CTRECiIR«CYUA, K) / IB 
CONTINUE 
TIME=SVALUE/3600.0 
СЕР=СОР+ Сир 
PRINT 360» CGPLCOP»CWTKP+TIME 
IF (CwPe6T.0.0) 320,310 
CONTINUE 
SVALUCZSVALUE*SVALPV/NDAT*10.0 
GO TO 5U 
CONTINUE 
PRINT 550 





330 
340 
320 


360 
370 


380 
390 
400 
410 
420 
430 
440 
450 


460 
470 
480 


490 


500 


510 
520 
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540 
290 


560 
510 
580 
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РЕТНТ 400% МТ, МР, МТ, МОАТ, ВТ. РВ. РОВ», ВАТТО» 2М.РМ 
PRINT SEO 

РАТНТ 410, РУОЁМХ.Н, ХМАХ », УМАХ, РУ 

PRINT 570 

PRINT 420» (QUID) Хи (Т1) s YW(CI) » Iz] 4 1 TOT) 

CALL ENoPLT 


FORMAT (4]1514F5.2+ F 7,0) 

FORMAT (Ғі0.0»5Х%Р 5.0е5Х.Ёрб.0) 

FORMAT (l6Xx++:"OBILIVY Cf FLUJL IN ZONE 1 = “fF10.3+/09 

1 l6X4eMOUBILITY CF FLUID IN ZONE 2 = #F)0.31/ 

2 lox+ MOBILITY CF ¿Ub IN ZONE 3 = *F10,.3s*/* 

3 16Xs+*MOBILITY CF FLUID IN ZONE 6 г 9Ғ10.3%/% 

4 l6Xs«40BILITY CF FLUID IN ZONE 5 = *F10.3) 

FORMAT (/16Х,610.2,8Х.Е10.2%8Х.610.2%125Х.Е10.2) 

FORMAT (7//1 6K a *GAS PRODUCED OUR RODUCED WATER PRODUC 
JED TIME IN HOURSt9/0]16X++CUBIC FEET RARRELS 
e BARRELS#) 

FORMAT ( /16X,F10e145X,F102135X*F5«1) 

FORMAT (4F 7.2517) 

FORMAT ( /16Х,415,10Е10,.2) 

FORMAT (/16Х%5Е10.3) 

FORMAT ( /16Х%610.3%5Х%Е10.3%5Х:Е10%3) 

FORMAT (3Ғ10.8) 

FORMAT (3F10«8) 

FORMAT (//16X*:17HSLOPE 11КЕ ONE 9 F9«7*3X9165HLINE CONST, = Fl 
10.8) ` 

FORMAT (16Х,еМІМІмУМ INJECTION RATE s *E10.2»/*16Y **MINIMUM PRO 


1DUCTION RATE = E1062 0/916Xe%wELL PATTERN AREA = Е] 0.2.1. 16Х 
2 DISTANCE BETWEEN INJECTION WELL AND PROVUCTION WFLL = #F10,2) 


FORMAT (16X,*^4MOUNT OF OIL IN ZONE 2 = #Е5,3:/ 16У AMOUNT CF O 
LIL IN ZONE 3 = ?F5,3,/:16Xs.eAMCUNT OF OIL IN ZONE 4 = 4F5,3) 
FORMAT (/16X** AVERAGE DENSITY OF RES, (RHO) = 9F10,4^*/516x«4*A'^QU 


INT OF OIL CONSUMED AS FUEL = 2Ғ10%49%/%16Х% <АУЕНАСЕ RES, OIL CONTEN 
eT GAL« PER TON z *F10«4) 

FORMAT (/16x»*RESERVOIR PORE VOLUME = %р9,2,%/ 916 УҰЗЭТКІСУМЕС6 0 
IF RESERVOIR = eF6,2./*16X,*MAXIMUM DIST IN A DIRECTION = еҒо,2,/% 
2 16x+*PAXÍMUM DIST IN Y DIRECTION = *F6.2+/9016X+*INJECTED PORE VO 
GLUES = #10) 

FORMAT (/16X,F NUMBER OF STREAM TUBES 7 *I54,/*16X**9 NUMBER OF PROD 
JUCTION WELLS = *1[5s/,.16Xe*NU'IRER OF INJECTION WELLS = “IE./416X+“N 
eUAT = 4]51/4016x9*RADIUS CF WELL BORE IN FEET = 9F5.24./*16X«.opOROSI 
ЗТҮ Ор РЕСЕНУОІН = #50 21/1 16%, HCOMPFRESSIBILITY FACTOD = °К5,2,,/,]6 
4X, *RESERVÜUIR PRESSURE IN ATMOSPHERES = %F700) 


FORMAT (2Ғ10.%) 

FORMAT (9F5.3) 

FORMAT (//16X*17H9LCPE LINE TWO 7 F9«7*3X»*14HLI^E CCNST, = F] 
16,9) 

FORMAT (4Ғ10.3) 

FORMAT (/18Х%%М5Т NP NI NCAT RI FMB POR RATI 
10 zM рме) 

FORMAT (//16X**PVOL MX H ХМАХ YMAX [РҮ AMB») 

FORMAT (//16X** WELL КАТЕ Xx LOC Y LAC +) 

FORMAT (3F 1004) 


END 
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APPENDIX C 
NOMENCLATURE 


Cross sectional area, L 


Constant 

Factor 

Thickness, L 
Permeability, L/ 
Length, L 

Moles, dimensionless 
Pressure, м/т? 
Injected volume, ү 
Injection rate, L°/T 

Universal gas constant (per mole), ML?/Tt 
Slope constant 

saturation, dimensionless 

Temperature, t 

Time, T 

Mux, L/T 

Volume, us 
Velocity, L/T 

Compressibility factor, LT4^/M 


Greek Letters 


D 
À 


Conductivity ratio, dimensionless 
Fluid mobility, L°T/M 
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u Viscosity, M/LT 

$ Potential, M/L?T? 

ф Porosity, dimensionless 

TT Constant = 3.1416 

p Density, M/L? 

Subscripts 

B Bulk 

b Burn 

eff Effective 

g Gas 

i Counter for number of zones 
К Counter for number of zones 
m Average 

n Maximum number of zones 
res Value at reservoir conditions 
S Standard fluid 

5С Value at standard conditions 
dE Total 

W Water 

X X-coordinate 

y y-coordinate 

2 Z-coordinate 

2 Zone 2 

3 Zone 3 

4 Zone 4 

9 Zone 9 





ADDRPO 


AMTOIL 
AREA 
BONE 
BTWO 
CO 

СОР 
CIR 


ery 
CW 
CWP 


CY 
DT 


FACB, 2, 3,4 


FACT2, FACT3 


FIVMB 
FLUXM 


FLUXMN 
FMB 
FORMB 
FUEL 
GALPTN 
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Total fluid accumulation for streamline under inves- 
tigation, L?/T 


Amount of oil in streamline under investigation, is 
Surface area bounded by production wells, 1,2 
Value of C in Equation 3.8 

Value of C between FLUXM and FLUXMN 
Cumulative oil produced by a streamline, L? 
Cumulative oil produced, field-wide, Lj? 


Pattern conductivity ratio for producing well, 
dimensionless 


Conductivity of individual flow tube, dimensionless 
Cumulative production of injected fluid 


Cumulative production of injected fluid field- wide 
basis, L 


Counter for the conductivity of individual flow tubes 
Time required to move one step on a streamline, T 


Denominator within the summation of the velocity 
equations for each zone, L4 


Term computed by Equations 3.37 and 3.41, 
dimensionless 


Mobility of fluid in zone 5, L°T/M 


Transition point of relation between air flux and 
advancement rate of burn, L?/T 


Minimum air flux required to move burn, L?/T 
Mobility of injected fluid, L?T/M 

Mobility of fluid in zone 4, L?T/M 

Amount of oil in place that is consumed by the burn 


дуегаре gallons of oil per ton of reservoir material, 
Г/М 





GASPRD 
GRATIO 
GRATI2, 3, 4 
H 

I 

IB 

IPV 

ITOT 


MULT 
NDAT 


NI 

NI1 

NP 

NSL 

NST 
NTOT 
РВ, 2, 3, 4 
PI 

PKST 

PM 

POR 

PP 
PVOLMX 
Q 


1 
Cumulative gas produced on a field-wide basis, L?/T 
Amount of gas in zone at production well 
Amount of gas in each zone 
Thickness of reservoir, L 
Counter for number of wells 
Counter for desired streamlines 
Number of pore volumes injected 
Number of image wells 
Step counter along streamline 
Streamline counter 
Counter for number of wells 


Number of even increments of SVALPV at which 
production history will be accumulated 


Number of injection wells 

First injection well 

Number of production wells 

Number of streamlines emanating from a well 
Number of streamtubes emanating from a well 
Total number of real wells 


Pressure of leading edge of each zone, M/LT? 


Pressure on well bore of injection well, м/т“ 


Constant preceding summation in Equation 3.5 
Average reservoir pressure, M/LT® 


Reservoir porosity, dimensionless 

Р 2 
Pressure on well bore of a production well, M/LT 
Reservoir pore volume 


Volume of fluid, L? 





QMIN 
QMINP 
QW 
RAD 


RATIO 
RATIO2, 3, 4 
RHO 

RI 

SCL 

SDTST 
SLOPEO 
SLOPET 


SMB 
SVALPV 
SVALUE 
SVELXB, 2,3,4 
SVELYB, 2,3,4 
TBT 

TIME 

TMB 

TSUM 

TSUMI 

VELB 
VELFAC 


VELSA 
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Minimum injection rate, L?/T 
Minimum production rate, [| т 
Volume of injected fluid, n 


Distance between point of interest and production 
well, L 


Portion of fluid in zone at production well 

Portion of fluid in zones 2,3,4 

Density of reservoir, M/L? 

Radius of well bore, L 

scaling factor to reduce scale of well 

Distance between production well and injection well 
Value of SL in Equation 3.8 


Value of SL when air flux is between FLUXM and 
FLUXMN 


Mobility of fluid in zone 2 

(Pattern area)*(H)*(POR) /QMIN, T 
(IPV)*(SVALPV) /NDAT, T 

Velocity of each zone in x-coordinate direction, L/T 
Velocity of each zone in y-coordinate direction, L/T 
Time of breakthrough of driving fluid, T 

Time since start of injection of recorded production, T 
Mobility of fluid in zone 3, L°T/M 

Cumulative time since start of injection, T 

Counter for cumulative time since injection, T 
Velocity of burn, L/T 


Ratio of velocity of burn to velocity of air across 
the burn 


Air flux, L/T 






VELSAR 
VELSB, 2, 3, 4 


VELTB, 2, 3, 4 


VELTMX 
VELXB, 2, 3, 4 
VELYB, 2, 3, 4 
VKST 


VOLO2, 3 
WATPRD 
WRATIO 
WRATI2, 3, 4 
X 

XEB, 2,3,4 


XFB, 2,3,4 


XMAX 
XW 

Y 

YEB, 2,3,4 


YFB, 2,3,4 


YMAX 
YW 
ZM 


Go 


Velocity of the air across the front of each zone along 
the streamline, L/T 


Velocity of the leading edge of each front along the 
streamline, L/T 


The maximum velocity of the above four, L/T 
Velocity of each zone along the x-coordinate, L/T 
Velocity of each zone along the y-coordinate, L/T 


Constapt preceding summation in Equation 3.6 and 
STALE 


Amount of oil in zone 2 and 3 
Cumulative water produced, E 
Amount of water in zone at production well 
Portion of zone 2, 3, 4 that is water 


X-coordinate of fluid interface along streamline, L 


X-coordinate of each zone along streamline L that 
is to be stored, L 


X-coordinate of each zone along streamline L to be 
used in computations, L 


X-coordinate of maximum point of interest, L 
X-coordinate of well locations, L 
Y-coordinate of streamline L on well bore, L 


Y-coordinate of each zone along streamline L to 
be stored, L 


Y-coordinate of each zone along streamline L to 
be used in computation, L 


Y-coordinate of maximum point of interest, L 
Y-coordinate of well locations, L 


Compressibility factor of driving fluid, LT%/M 
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12. 
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